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We introduce a class of Markovian quantum master equations, able to describe the dissipative dy¬ 
namics of a quantum system weakly coupled to one or several heat baths. The dissipative structure 
is driven by an entropic operator, the so called modular Hamiltonian, which makes it nonlinear. 

The generated Modular Dynamical Semigroup (MDS) is not, in general, a Quantum Dynamical 
Semigroup (QDS), whose dynamics is of the popular Lindblad type. The MDS has a robust ther¬ 
modynamic structure, which guarantees for the positivity of the time evolved state, the correct 
steady state properties, the positivity of the entropy production, a positive Onsager matrix and 
Onsager symmetry relations (arising from Green-Kubo formulas). We show that the celebrated 
Davies generator, obtained through the Born and the secular approximations, generates a MDS. 

By unravelling the modular structure of the former, we provide a different and genuinely nonlinear 
MDS, not of QDS type, which is free from the severe spectral restrictions of the Davies generator, 
while still being supported by a weak coupling limit argument. With respect to the latter, the 
present work is a substantial extension of fli- 
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Quantum Markovian master equations 3j-|f| are be¬ 
ing used since half a century to model the dynamics 
of a quantum system in contact with a particle reser¬ 
voir 0 , 0 , 0 . If on the one side they are justified on the 
statistical mechanical grounds of an underlying purely 
Hamiltonian system in a scaling limit 00 , their interest 
often stems from their own structural properties, partic¬ 
ularly with resp ect to thermodynamic requirements, on 
the other l9Hlll|. A famous example is the Davies master 
equation [3|, obtained through the Born plus secular ap¬ 
proximations Q and giving rise to a Quantum Dynamical 
Semi grou p through a generator of the popular Lindblad 
kind [12j. As many other known generators, the Davies 
generator can (i) describe the Hamiltonian evolution of a 
small subsystem interacting with a particle reservoir at 
equilibrium in the weak-coupling limit 0,[l3| . The Davies 
generator turned out to be the most successful and pop¬ 
ular choice, as it comes with the additional advantages 
of guaranteeing (ii) the complete positivity of the time 
evolved state [3, (hi) the correct steady state 0, (iv) 
the positivity of the entropy production, (v) a positive 
Onsager matrix fl5j and (vi) Onsager symmetry relations 
for heat baths in the linear regime (arising from Green- 
Kubo formulas) 0EM3- That all these properties are 
fulfilled for finite values of the coupling to the bath, and 
not just asymptotically in the scaling limit, has given the 
status of a physically robust quantum dissipative dynam¬ 
ics per se to the Davies generator. 

However, Davies himself in his foundation papers 00 
pointed out that his generator suffered from severe spec¬ 
tral limitations, and could not even be defined for un¬ 
bounded subsystems. In [30 this problem was stud¬ 
ied, and a cut-off dependent Lindblad generator was pro¬ 
posed. The generator was supported by the same weak- 
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coupling argument as for the Davies case, and unlike the 
latter it did not suffer from any spectral restriction. How¬ 
ever, nothing could be said in general for properties (iii)- 
(vi) above, at finite coupling to the bath. Therefore, 
a master equation able to extend the Davies generator 
beyond its spectral limitations, while preserving (i)-(vi) 
above, is still lacking. Moreover, even for bounded sub¬ 
systems, the important question remains, as to weather 
the Davies master equation is unique: Is there really no 
other possible Markovian dynamics left, which guaran¬ 
tees (i)-(vi) above? 

In this paper, we shall present a class of nonlinear 
Markovian master equations, describing the dissipative 
time evolution of out of equilibrium excitations of a quan¬ 
tum system coupled to one or several heat baths. The 
scattering event governing the irreversible dynamics is 
driven by an entropic operator at the origin of the non¬ 
linearity, known as the modular hamiltonian. For this 
reason we say that our dynamics generates a Modular 
Dynamical Semigroup (MDS). We prove a thermody¬ 
namic analogue of complete positivity (the latter notion 
is not even available for a nonlinear dynamics), namely 
that (ii’) the relative entropy operator is bounded at all 
positive times, even when an ancilla system is included. 
Properties (ii’-vi) are all consequences of the modular 
structure, as we show. MDSs are at variance with the 
linear QDSs. We find that the only QDS which is also 
a MDS is that given by the Davies generator. Unrav¬ 
elling the modular structure of the latter immediately 
leads to a possible nonlinear MDS extension, able to de¬ 
scribe unbounded systems and systems in the thermo¬ 
dynamic limit. That the generalization is not unique 
gives substantial freedom for engineering novel dissipa¬ 
tive quantum phenomena through a simple, Markovian 
master equation. Being now able to properly describe 
systems in the weak coupling regime, this work is a struc¬ 
tural extension of 00. As therein, more general envi¬ 
ronments can be taken into account, although here we 
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only consider heat baths for sake of exposition. 

We take H = 1 and let the states of the quantum system 
be described by strictly positive and normalized trace 
class operators p over a Hilbert space H. Observables 
A are bounded operators on 'H. The system is put in 
contact with a heat bath at inverse temperature 0 and 
otherwise undergoes a free hamiltonian dynamics gener¬ 
ated by a Hamiltonian Hg (potentially containing renor¬ 
malization terms). The irreversible influence of the heat 
bath on the system is implemented through a bounded 
scattering operator Q x depending on A S [0,1] and on 0. 

We let AS = — In p — 0Hg be the relative entropy op¬ 
erator. We define a (0, p)-dependent modular dissipative 
bracket [[•, •]]/?, P , as the positive bilinear form over the 
observables 

[[A,B]] PiP =£(cj? WVt]))# (i[Q x ,B])d\, 

. (1) 

constructed with the help of the hamiltonian evolution 

(‘(H) = e “(-/3- 1 inp)H e - it( - /3_lln ‘ ,) , (2) 

generated by the modular hamiltonian —S 1 lnp at the 
imaginary time t = i\0/2 (see e.g. [191120|). 

Through such brackets, and denoting p(A ) = tr(pH) 
when there is no room to confusion, we define the MDS 
Wt,t' to be the evolution map W t .t' ijh') = Pt associated 
to the master equation 

9tp{-) = P (i[H s , •]) + p ([[-, AS]]^). (3) 

We note that this is a generalization of the nonlinear 
quantum master equation presented in 00 , which is 
recovered in the particular case of Q x = Q not depend¬ 
ing upon A (with the understanding that extension of 
the present formalism to more general environments than 
heat baths can be made by letting the environment brack¬ 
ets in [l| depend on A). This is an important conceptual 
development, motivated by a much wider applicability 
regime of the resulting theory, as we shall see. It allows 
the bath to have a more general kind of influence the 
system, thanks to what could be named colored modular 
noise. Indeed, the evolution © at the microscopic imagi¬ 
nary time i\0/2, during the scattering event, is weighted 
by the coupling operator Q x . On the top of this colored 
modular noise one could accommodate for colored noise 
in the usual sense by extending © to many coupling op¬ 
erators Q x for a € I in some energy level set (as we shall 
later do). The interpretation of © is clearly that the 
relative entropy operator AS acts as a thermodynamic 
force driving the system to equilibrium, just as in 0 . 

For the well posedness of © it is sufficient to as¬ 
sume, as we do, that Q A ^ = Q 1 ^ x (which guarantees 
self-adjointness), and that [Q x ,Hg] is bounded. We fur¬ 
ther make the ergodicity assumption that Hg and Q x 
do not share common eigenspaces 0 m. This guaran¬ 
tees that all scattering channels are coupled, leading to a 
unique steady state, when the latter exists (in which case 
it is the equilibrium state pp, as seen from AS = cl). 


Note that the modular evolution (‘in © is still well 
defined in the thermodynamical limit [19| , where the pos¬ 
itive normalized functional p(-) implementing the state of 
the system can no longer be represented by a density ma¬ 
trix [20j]. As a direct consequence the MDS master equa¬ 
tion © easily extends to the thermodynamical limit, at 
least as long as AS is bounded throughout the evolution. 
Note that if H is finite dimensional, as we shall from now 
on assume to avoid technicalities, this is the case for just 
any strictly positive p > 0 . 

This leads us to the following positivity argument. 
First, © preserves normalization through <9 t p( 1) = 0. 
The only way for A S' to diverge is then for some eigen¬ 
value p of p to go to 0. When this happens, the en¬ 
tropy S(p) = p(—]np) decreases. But from d t S(p) = 
p([[— lnp, AS]]/ 3 jP ) and from [Q a ,AS] = [Q A ,lnp] +0(1) 
(due to the supposed boundedness of [Q x ,Hg]), we find 
that in this limit, the Von Neumann entropy rate is 
asymptotic to 

°{p) = P ([[AS, AS]]^) > 0, (4) 

which is positive (because of the positivity of ©, thus 
guaranteeing that AS stays bounded throughout the evo¬ 
lution. 

We wish at this point to discuss the positivity property 
of the MDS in comparison to the notion of complete pos¬ 
itivity BUI that cannot even be defined for our nonlin¬ 
ear dynamics (which is not even of mean-field type [ 2 ^ 1 . 
We take the physical essence of the disputed notion (see 
e.g. t 23j and references therein) to be that the presence 
of a external ancilla quantum system S', non-interacting 
with S, cannot spoil the positivity of the evolution. From 
our perspective, the most natural implementation of such 
concept is that the total hamiltonian and the total cou¬ 
pling operator of S + S' should not couple degrees of 
freedom of S and S'. We therefore take them of the 
form Hg ® lgi + Is <g> Hgf and Qg <g> I 5 / + Is ® Qg, re¬ 
spectively. Then the MDS of S naturally extends to an 
MDS of S' + S', which is again positive (and so it stays 
in the limit of vanishing Hg> and Qg,, provided that er¬ 
godicity is guaranteed throughout the limit). This is why 
we can safely state that MDS are completely positive in 
this sense. We believe that our alternative and inequiv¬ 
alent notion of complete positivity, although trivial and 
innocent, can shine light on the problematic conventional 
definition, which itself exhibits such surprisingly drastic 
physical consequences on the system S [A| . 

Actually there is a deeper thermodynamical meaning 
to © : by defining the mean entropy flux from the system 
to the reservoir as Jg{p) = —0d t p{Hg), one sees that © 
gives nothing else than the entropy production a(p) = 
d t S{p) + Jg{p) = dtp(AS), analogously to what already 
remarked in |2j. Since er(p) is strictly positive unless 
AS = cl for a c-number (due to our ergodic hypothesis), 
this implies both the uniqueness and the thermodynamic 
stability of the steady state p + = pp. 

In the belief that the ergodicity assumption should im¬ 
ply that the steady state is unique, when it exists, also 
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for the case of N heat baths at inverse temperatures Pj, 
j = 1 ...,7V, we now define the Onsager matrix (l5 j in 
the linear regime, where the entropy fluxes are linearly 
proportional to the thermodynamic forces Xj = pj — p 
(for a reference inverse temperature P). We prove that 
the Onsager matrix is positive and give a Green-Kubo 
formula jl6l . 0 for its coefficients (thereby obtaining 
the Onsager symmetry relations). These ideas and tech¬ 
niques combine the analyses in [24l . Erl for a time-driven 
setting, and in 0 for the Davies case. 

To this end, we consider the MDS d t p = D(p) = 
p(i[Hs, •]) + Djpjip) where each dissipative con¬ 
tribution Djp.(p) = p ([[■, A5j]]^ jp ) is given in ([3]) with 

scattering operators Qj.g., and relative entropy operator 
A Sj = A= In p — In pp r Note that pp is the only 
steady state at equilibrium pi = ... = /3jv = P (as can be 
seen from the above entropy production argument) and 
that, under our ergodic assumption, p«. is the unique 
steady state of D d g, (•). 

Using K p A = fij p x Ap 1 ~ x dX (which admits an in¬ 
verse) we compute the linearization D d g{p) = D d gp + 
o(p — pp) at the equilibrium [26| to be 

Dj,p P = -PP([ h K ppP]]p,p p )- ( 5 ) 

We see that D d g is purely dissipative, e.g. D d gp = ujp 
only has solutions with w < 0, and that pp is the 
unique steady state of D d g . Indeed by tracing the eigen¬ 
value equation against p, the right hand side reads 
ujtv{pK~ d p), and the trace is positive, as follows from 
a calculation in the H$ basis. The left hand side is the 
negative of 

Pp([[K^p,K^p% t J (6) 

and is thus strictly negative unless p = pp. The argument 
follows. Note that D d g need not be of Lindblad type. 

The entropy production at the steady state is cr(p+) = 
(J j{p+) = ]CjLi XjJj where the mean fluxes Jj = 
p + ([[H s ,AS k (p + )]]^ tp+ ) = Ef=i X k L jk (J3) + o(X) de¬ 
fine the Onsager coefficients Lj k (P). From the above lin¬ 
earization we deduce that the entropy production cr(p+) 
next to equilibrium is a bilinear form in p + — pp, whose 
terms <Jj{p + ) are precisely given by (O with [3 = pj. In 
particular pp is a local minimum for cr(p+), which implies 
that the Onsager matrix is positive definite (over the re¬ 
als). To prove the Onsager relations L k j{P) = Lj k {P) we 
need a notion of detailed balance next to equilibrium. To 
achieve this, we note that we can safely trade the mean 
fluxes Jj with their linearization Jj made at pj (according 
to ©), and thus limit ourselves to analyze the detailed 
balance properties of D d = D d p- a t equilibrium. 

Defining D = —i[Hs, ■} + D d and D* on BpH) by duality, 
we note that at equilibrium 

D*=K~ d DK Pf> . (7) 


We remark that this detailed balance property is con¬ 
tained in the modular structure of the MDS, and does not 
come a priori from the KMS property fijjj of some bath 
correlation function Q. Of course, these two notions 
can be linked through our specific choice of the coupling 
operator Q x , as we shall see, but the conceptual sepa¬ 
ration is evident. We are now in position to go through 
Lemma 1 of jp . with the only difference that we now get 
Dip{H s pp) = D d g(K Pfj H s ) = K pp Df *(H s ) at equi¬ 
librium. By defining the flux operators Jj = D^’g(Hs), 
and using the fact that D is dissipative (as can be sim¬ 
ilarly proven for D* as well) the Green-Kubo formula is 
obtained 

/•OO 

Ljk(P) = / ; 3k)pdt. (8) 

Jo 

The r.h.s. is written with the help of the Kubo scalar 
product ( A ; B)p = t r(A^K pp B) and the Heisenberg evo¬ 
lution Apt) = e D t A. The detailed balance property 0 
implies that JJ j(t) = K~ d e Dt {K pfj Jj), from which it fol¬ 
lows that (JJj(Z) ; lfc )/3 = (JJfe(Z) ; JJj)/ 3 , providing the On¬ 
sager symmetry relations. Note that the currents Jj need 
not commute with pp, so our formula is in general differ¬ 
ent from the one provided in [§], and is exactly the well 
known Kubo formula |T(1 . 

The analysis above puts the MDS on solid thermody¬ 
namical grounds. We shall now show that the set of inter¬ 
esting MDS is far from empty. The first notable example 
is the celebrated Davies generator Q. To see that it be¬ 
longs to our MDS class, take eigenoperators A v defined 
by [A v ,Hs\ = vA u , a bath spectral function h{v) > 0 
with the KMS property h{y) = e^"h(—v) 0, and define 

Q x = e~ x ^' J l' 1 \jh(v)A v . Then d3j) is generalized to the 
sum over all v of the corresponding brackets [[•, •]]£ , and 
is guaranteed to be self-adjoint by the natural generaliza¬ 
tion Qpp = to many coupling operators Q x . The 

resulting master equation reads 

dtp = -i[H s , p\+Y^ H") P) + A v pA]\ . 

V ' ' 

(9) 

This can be seen from J 0 ' e~ Xf3u p x [A„, F]p 1 ~ x dX = A u p— 
e~^ v pA v , as can be computed in the p basis, and then 
using the KMS property of h. The r.h.s. of Eq. m gives 
exactly the Davies generator fsj], describing the weak 
coupling limit of a subsystem S coupled to a thermal 
bath B at inverse temperature p through a hamiltonian 
H = Hs + + Hb , where h(y) is Fourier 

related to h(t) = {e lHBt <&e~ %HBt §)p and R = J2 V A V . 
In passing, we remark that the Davies generator arising 
from a more general coupling hamiltonian ^ R a ® $ a 
can just as well be cast into MDS form, through a trivial 
extension of the above with eigenoperators A a relative to 
frequencies v a . We note that the KMS property of the 
bath correlation functions is here embedded as a neces¬ 
sary condition to have the MDS well defined: This is the 
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natural generalization = Q _J' of the self-adjointness 
condition Q A ' f = Q 1 ~ x to many coupling operators Q x . 
That the linearization of the Davies generator according 
to © is again the Davies generator was already proven 
in [25j- Since in this case [1 j,Hs] = 0, our expressions 
for the Onsager coefficients agrees with those computed 
in 

The above unravels the thermodynamic and modu¬ 
lar structure of the Davies generator, and automatically 
gives different proofs and more symmetric expressions for 
the entropy production, the fluxes and the Onsager re¬ 
lations. We make the following remarks. The bounded¬ 
ness of the relative entropy operator is more relevant to 
a thermodynamic setting than complete positivity. The 
former is guaranteed by the modular structure in (f3j). 
while the latter by the Lindblad structure [l3| . Second, 
from the modular perspective, the linearity of the Davies 
generator, together with its complete positivity, appears 
to be a fortunate, accidental coincidence, which is due 
to the peculiar commutation relations of the eigenoper- 
ators A„ appearing in the coupling operator Q x . The 
eigenoperators exist however only when H$ has a purely 
discrete spectrum, which severely limits applications to 
systems without internal structure, namely quantum dots 
and harmonic oscillators. Moreover even for those, other 
alternative MDS could be of interest. 

These remarks led us to a previous work by one of 
us [l3j, where the eigenoperators A„ were traded with 
different scattering operators 

/ +oo 

e il/t ^/5(t,T) e iHst Re iHst dt, (10) 

-OO 

where v G R, S(t,T) is the normalized Gaussian at t, 
whose standard deviation T physically represents a col¬ 
lision time (which can depend on a coupling constant, in 
a weak coupling limit spirit). Contrary to the eigenop¬ 
erators, the A„’s are always well defined and bounded 
(as long as 0 < T < oo), independently of the spectral 
properties of Hs . Proceeding as before, we now take our 

coupling operators to be Q x = e -A/3l/ / 2 \Jh{v)A„, and 

with those, we define the dissipative brackets [[•, m ]\p tP like 
in m■ The resulting MDS master equation 

dtp(-) = p ( i[H s , •]) + P [[-, A S]] V p if> dvj (11) 

is still well defined (Aj) = A_ u implies Q At = Q^ A ) 
and does not suffer from any spectral restriction on the 
system Hamiltonian. However, the above MDS master 
equation will in general be genuinely nonlinear. It can 
be seen with analogous calculations as in 0 that for 
Hs with discrete spectrum it boils down to the Davies 
generator when T oo. As before, the analysis easily 
accommodates for more general coupling hamiltonians 
J2a ® < E’ a . It follows that (fill) is fully justified on the 
statistical mechanical basis of the weak coupling limit, 
and provides a genuine and thermodynamically robust 


MDS extension of the Davies master equation, able to 
describe even unbounded systems, possibly in the ther¬ 
modynamic limit. 

We also remark that when the spectral function h{v) 
is singular at v = 0 (white noise), then the only coupling 
operator involved in m is (possibly some n-th order 
derivative of) Q x at v = 0, which will not depend upon 
A. In this case, m gives back the nonlinear master 
equation in [l, 2j. This raises the observation that m 
should be regarded as a (nonlinear) quantum counterpart 
of the classical linear Boltzmann equation, taking the 
singular limit of the scattering rates of which (28j gives 
the (nonlinear) quantum counterpart @] of the classical 
Fokker-Planck equation. 

We finally want to address the problem of unicity of the 
MDS fulfilling the properties (i-vi) stated in the introduc¬ 
tion. Although we have already proven the existence of 
two MDS master equations for bounded systems (© and 
CD), many more possibilities remain. As an additional 
example, one could consider a single coupling operator 

Q x = e~ X ^ v ^ 2 \Jh(v)A„, and the resulting nonlin¬ 

ear evolution will still match the corresponding coupled 
hamiltonian evolution in the weak coupling limit. To see 
that, it is sufficient to take the time average map of the 
MDS master equation |3| and realize that the result is 
just © (a solid argument as to why this should imply 
agreement in the weak coupling limit can be found in [7j). 
However, while these MDS all agree in a limit regime, 
they will in general give (possibly qualitative) different 
predictions at small but finite values of the coupling to 
the bath, which is where master equations are used in 
applications. 

To summarize, we have proposed the class of MDS 
master equations, describing the Markovian, nonlinear 
dissipative dynamics of a quantum system interacting 
with one or more particle reservoirs at equilibrium. The 
MDS is structurally and thermodynamically robust, as 
opposed to a general QDS. From this perspective, we 
find that the Davies master equation generates the only 
MDS which, accidentally, also happens to be a QDS. This 
is just due to the peculiar character of the correspond¬ 
ing coupling operators. We have used the MDS struc¬ 
ture of the Davies generator to extend the latter to a 
more general, nonlinear MDS master equation, which is 
free form the severe restrictions of the former and can 
even describe unbounded systems, possibly in the ther¬ 
modynamical limit. Other possible MDSs agree with the 
Davies MDS in the weak coupling limit, while differing 
from it at finite values of the coupling to the bath. The 
origin of the nonlinearity of the MDS lies in the intrinsic 
entropic dependence of the scattering event, and is very 
different in nature from any mean field approach, where 
nonlinearities are normally expected to arise. The MDS 
master equation relative to state-independent coupling 
operators should be regarded as the quantum counter¬ 
part of the classical linear Boltzmann equation. 
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